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Abstract 

We investigate the prospects for the discovery at the CERN Large Hadron Collider or at the 
Fermilab Tevatron of neutral Higgs bosons through the channel where the Higgs are produced 
together with a single bottom quark and the Higgs decays into a pair of tau leptons, bg — > b<fP — > 
6t + t~, (fp = h°, H°, A . We work within the framework of the minimal super symmetric model. The 
dominant physics background from the production of 6r + r~, jr + r~ (j = g,u,d, s,c), bbW + W~ , 
W + 2j and Wbj is calculated with realistic acceptance cuts and efficiencies. Promising results are 
found for the CP-odd pseudoscalar (^4°) and the heavier CP-even scalar (H°) Higgs bosons with 
masses up to one TeV. 

PACS numbers: 14.80.Cp, 14.80.Ly, 12.60.Jv, 13.85Qk 
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I. INTRODUCTION 



In the minimal supersymmetric standard model (MSSM) the Higgs sector has two 
doublets, <pi and 0o, which couple to fermions with weak isospin t 3 = — 1/2 and £3 = 
+ 1/2 respectively |2|. After spontaneous symmetry breaking, there remain five physical 
Higgs bosons: a pair of singly charged Higgs bosons H ^ , two neutral CP-even scalars H° 
(heavier) and h° (lighter), and a neutral CP-odd pseudoscalar A . The Higgs potential is 
constrained by supersymmetry such that all tree-level Higgs boson masses and couplings are 
determined by just two independent parameters, commonly chosen to be the mass of the 
CP-odd pseudoscalar (Ma) and the ratio of vacuum expectation values of the neutral Higgs 
fields (tan/3 = V2/V1). 

At the CERN Large Hadron Collider (LHC), gluon fusion (gg ^0,0 = h°,H°, or A ) 
is the major source of neutral Higgs bosons in the MSSM for tan/3 less than about 5. If 
tan/3 is larger than 7, neutral Higgs bosons are dominantly produced from bottom quark 
fusion bb — > jD, 0, H, H, @] • Since the Yukawa couplings of (pbb are enhanced by 1/ cos /3, the 
production rate of neutral Higgs bosons, especially the A or the H°, is enhanced at large 
tan p. 

Recently, it has been suggested that the search for a Higgs boson produced along with a 
single bottom quark with lar ge t ransverse momentum (p?) , where the leading order subpro- 
cess is bg — > b<f) jl, 0, 10, 11, 12], could be more promising than the production of a Higgs 



boson associated with two high p T bottom quarks [IOJ where the leading order subprocess 



is gg — > bb<f> [3|, [l3|, [14|, [l5|, [l6j. This has already been demonstrated to be the case for 



the u + u decay mode of the Higgs bosons [17|. For a large value of tan/3, the r + r decay 



mode [18|, [l9| is also a promising discovery channel for the A and the H° in the MSSM 
because the branching fraction for Higgs decay into tau leptons is greater by a factor of 
(m T /m At ) 2 ~ 286. The downside is that unlike muons, tau leptons can only be observed 
indirectly via their hadronic or leptonic decay products. 

In this article, we present the prospects of discovering the MSSM neutral Higgs bosons 
produced with a bottom quark via Higgs decay into tau pairs. We calculate the Higgs signal 
and the dominant Standard Model (SM) backgrounds with realistic cuts and efficiencies and 
evaluate the 5<r discovery contour in the (Ma, tan /3) plane for the LHC and for the Tevatron. 



II. THE PRODUCTION CROSS SECTIONS AND BRANCHING FRACTIONS 

We calculate the cross section at the LHC for pp — > b<p + X and at the Tevatron for pp — > 
b<t> + X (<f) = H°,h°,A°) via bg — > b<p with the parton distribution functions of CTEQ6L1 
[20( | . The factorization scale is chosen to be jfj], |2lj]. In this article, unless explicitly 

specified, b represents a bottom quark (b) or a bottom anti-quark (b). The bottom quark 
mass in the cf)bb Yukawa coupling is chosen to be the next-to-leading order (NLO) running 
mass at the renormalization scale /ar, mb(fiR) |22| . and it is calculated with m&(pole) = 4.7 
GeV and NLO evolution of the strong coupling [23(. We have also taken the renormalization 
scale for the production processes to be M^/4, which effectively reproduces the effects of 
next-to-leading order [Io|. Therefore, we take the K factor to be one for the Higgs signal. 

The cross section for pp — > b<f) — ► br + r~ + X can be thought of as the Higgs production 
cross section a(pp — > b<p + X) multiplied by the branching fraction of the Higgs decay into 
tau pairs B(<p — > t + t~). When the bb mode dominates Higgs decays, the branching fraction 
of — > t + t~ is about / ("Sm^M^) + m 2 ) where m^M^), the running mass at the scale 
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M^, is used in the decay rates. This results in a branching fraction for A — > t + t of ~ 0.1 
for M A = 100 GeV. Thus for tan/5 > 10 and M A > 125 GeV, the cross section of bA° or 
that of bH° is enhanced by approximately tan 2 (3 and the branching fraction of Higgs decay 
to tau pair is close to 10%. 



III. TAU DECAY AND IDENTIFICATION 

Tau leptons can decay either purely leptonically, t~ — > £~i?£U T , with a branching ratio 
of around 18% for each mode I = e, fj,, or they can decay into low-multiplicity hadronic 
states and a v T with a branching ratio ~ 64% [24]]. Therefore, for a r + r~ pair, the most 
likely scenario is one decaying leptonically and the other hadronically, which has a combined 
branching ratio of 46%. Also, the presence of an isolated lepton in the final state is useful 
in triggering the event and reducing backgrounds. Hence, we use this "lepton + r-jet" 
signature in our study. 

We model hadronic tau decays as the sum of two-body decays into 7rz/ T , pv T and a\v T with 
branching ratios given in the literature [24| . The tau is assumed to be energetic enough that 
all its decay products emerge in approximately the same direction as the tau itself. This 
manifests itself in the so-called "collinear approximation" which we use for both leptonic 
and hadronic decays. The approximation is confirmed to be accurate by comparison with 
an exact matrix element simulation for tau decay. 

In Figure 1, we present the transverse momentum distribution [dajdpT] for the bottom 
quark (b), or the lepton (I) or the tau hadron (j T ) from tau decays, for the Higgs signal 
pp — > bA° — > br + r~ —> Mj T + X. In addition, we show the px distribution for b, £, or j T 
from the SM background bg — > br + r~ (Drell-Yan). We have required pr{b) > 10 GeV and 
\r]b\ < 2.5. The purpose of this figure is to show these cross sections before any other cuts 
have been applied. 



The ATLAS collaboration has studied identification efficiencies of r-jets in detail [25 
Based on this we use an overall efficiency of 26% over 1- and 3-prong decays with a corre- 
sponding cut, px(h) > 40 GeV for the hadron h = n, p, a±. This also corresponds to a mistag 
efficiency of 1/400 for non-r (i.e. QCD) jets. Rejection of jets from b quarks is higher, with 
only 1 in 700 being mistagged as rs. The transverse momentum cut on the lepton from tau 
decay is weaker, with pr{£) > 20 GeV. Both the hadron and lepton are required to be in 
the central rapidity region \r]\ < 2.5. The acceptance cuts as well as tagging and mistagging 
efficiencies for the Fermilab Tevatron will be discussed in Section VII. 



IV. HIGGS MASS RECONSTRUCTION 

The Higgs mass can be reconstructed indirectly, using the collinear approximation for r 
decay products and the missing transverse momentum 2-vector, ^ T . Taking xg,Xh to be 
the energy fractions carried away from the decays by the lepton and hadron respectively, we 
have: 

(- - 1)P" T + (- - 1)Pt = # T (1) 

Xt Xh 

This yields two equations for xi and Xh which can be solved to reconstruct the two original 
r 4-momenta p^ = p$/xi, Ph/ X h- Thus M 2 = (pi/xi +Ph/xh) 2 - Physically we must have 
< xe, Xh < 1, and this provides a further cut to reduce the background. 
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(b) b g — > b x + t" 
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The transverse-momentum distribution for (a) the Higgs signal from bg — > bA° — > 
Mj + X with Ma = 200 GeV and tan (3 = 10 as well as for (b) the physics background 
-> br~r + — > Mj + X,j = n,p, or a±. In the three curves, px refers to the transverse 



momentum of the b— quark or the lepton or the tau-jet. 



Measurement errors in lepton and r-jet momenta as well as missing transverse momentum 
give rise to a spread in the reconstructed mass about the true value. Based on the ATLAS 
and the CDF specifications we model these effects by Gaussian smearing of momenta: 



AE 0.50 



E 



0.03 



for jets (with individual terms added in quadrature) and 

AE 0.25 



E 

AE 
~E 



'E 
0.15 



E 



0.01 (LHC) 
0.01 (Tevatron) 



(2) 

(3) 
(4) 



for charged leptons. 

We find that in more than 95% of the cases, the reconstructed mass lies within 15% of 
the actual mass. Therefore we apply a mass cut, requiring the reconstructed mass to lie in 
the mass window M^± AM TT , where AM TT = 0.15M^ for an integrated luminosity (L) of 30 
fb _1 and AM TT = O.2OM0 for L = 300 fb _1 . This cut is actually rather conservative because 
for larger Higgs masses, more than 90% of the reconstructed masses are within 5 — 10% of 
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M^. We note that improvements in the discovery potential will be possible by narrowing 
AM TT if the r pair mass resolution can be improved. 

Figure 2 shows the invariant mass distribution of the tau pair for the Higgs signal pp — > 
bA° — > br + r~ + X via bg — > bA°, and for the tau pair from the SM Drell-Yan background 
bg — > br + r~ . We have calculated the Higgs signal in two ways: (a) with the narrow width 
approximation 

a(pp -> bA° -> br + r- + X) = a(pp -> bA° + X) x -> 6r + r-) 

and (b) the full calculation er(pp — > — > br + r~ + X) with a Breit-Wigner resonance via 
bg — > — > br + r~. In this figure we have applied all acceptance cuts discussed in the 
next two sections except the requirement on invariant mass. We note that with energy- 
momentum smearing, the cross section in the narrow width approximation agrees very well 
with that evaluated for a Breit-Wigner resonance. 
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FIG. 2: The invariant-mass distribution distribution, da/dM TT (pp — > br + T~ — > b£j T + + X), 
for the Higgs signal from bg — > bA° — > t~t + — > Mj T + X with Ma = 150 GeV or 200 GeV and 
tan/3 = 10 as well as for the physics background from bg — > br~r + — > b£j T + X,j T = n,p, or a\ 
(dash). We calculate the Higgs signal in the narrow width approximation (short dash) and with a 
Breit-Wigner resonance from a Higgs propagator (dash-dot). 
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V. THE PHYSICS BACKGROUND 



From the above discussion, the signal we are looking for is 6-jet (b) + lepton (£) + r-jet (j) 
+ Ifirp + X, where Ifi T = missing transverse energy ~ $ T = missing transverse momentum. 
The dominant physics backgrounds to this final state come from: 

(i) Drell-Yan processes: pp — ■> jZ*/j*+X — > jr + r~+X, j = u, d, s, c, b, g. Approximately 
60 — 70% of the DY contribution arises from the subprocess bg — > br + r~ . 

(ii) Top Production (gg,qq — > tt — > bbW + W~): This can contribute in several ways 
depending on how the Ws decay. In order of highest to lowest importance, the relevant 
channels are the following, (a) brvblv: One W decays into tv t with the r decaying 
hadronically while the other W provides tv^. (b) brvbrv: We can have both Ws 
decaying into tv t with one tau decaying leptonically and the other hadronically. (c) 
bivbjj: In this case we can have one W decay leptonically while the other W decays 
into jets (W — > qq'). We now have four possible jets in the final state i.e., 2 6's and 2j 
and one of them is tagged as a b quark while one of the other is mistagged as a r-jet. 
(d) btvbiv: We can have both Ws decay leptonically. Then we have two b quarks 
in the final state, one of them is tagged as a 6 — jet while the others is mistagged as 
a r-jet. (e) brubjj: Finally we can have one W decay into tv t with the r decaying 
leptonically while the other W decays into jets (W — > qq'). We now have four possible 
jets in the final state i.e., 2 6's and 2j and one of them is tagged as a b quark while 
one of the others is mistagged as a r-jet. 

(iii) tW Production (bg — > tW — > bW + W~): This is very similar to the top quark pair 
production just discussed. In order of decreasing importance, the relevant channels 
are as follows, (a) bruiu: One W decays into rv T with the r decaying hadronically 
while the other W provides Iv^. (b) brvrv: We can have both Ws decaying into tu t 
with one tau decaying leptonically and the other hadronically. (c) btvjj: In this case 
we can have one W decay leptonically while the other W decays into jets (W — > qq'). 
We now have one b and 2j with one of the light quarks mistagged as a r-jet. (d) brujj: 
Lastly we can have one W decay into tv t with the r decaying leptonically while the 
other W decays into jets (W — > qq'). Again, we have one b and 2j with one of the 
light quarks mistagged as a r-jet. 

(iv) W + 2j processes: pp —> W + 2j + X with the subsequent decays W — > iui, £ = e, fi 
or W —> tu t with the r decaying leptonically. Here, one jet is tagged or mistagged as 
a b quark and the other mistagged as a r-jet. 

Due to the huge cross-section for pp qqg with q = b,c, it is also pertinent to check 
that heavy quark semi-leptonic decays such as b — > civ do not overwhelm the signal. We 
find that this background is effectively cut to less than 10% of the dominant background at 
all times by an isolation cut on the lepton \rj(£, j)\ > 0.3, the large rejection factor for non-r 
jets, and the requirement Ifi T > 20 GeV. 

For the lower integrated luminosity (L) of 30 fb _1 , we require pr{b,j) > 15 GeV and 
\r)(b,j)\ < 2.5. The 6-tagging efficiency (e&) is taken to be 60%, the probability that a c-jet 
is mistagged as a 6-jet (e c ) is 10% and the probability that any other jet is mistagged as a 
6-jet (ej) is taken to be 1%. For the higher luminosity L = 300 fb" 1 , we take e& = 50% and 
p T (b,j) > 30 GeV @. 
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In order to improve the signal significance we also apply a cut on the transverse mass 
of tut^i $ T ) < 30 GeV. Using the definition of transverse mass given in 27j we find that 
this is very effective in controlling the W + 2j and ti backgrounds. In addition we require 



25, 26 , for the 



, r — jet) < 170°, as suggested by ATLAS and CMS collaborations 
reconstruction of the Higgs mass as the invariant mass of tau pairs. 

We have applied a K factor of 1.3 for the DY background [28J1, a K factor of 2 for ti 
2§|, a K factor of 1.5 for tW 



for bq — > Wbq, q = u.d, s.c [32 



30 



a K factor of 0.9 for W + 2j [31] , and a K factor of 2 
to include NLO effects. In order to further cut down the 
ti background, we apply a veto on events with more than 2 jets in addition to the b and r 
jets. This is very effective because, in ti + X production, nearly 50% of events have at least 



one gluon from initial or final state radiation that passes px > 15 GeV and \rj\ < 2.5 [291 ] . 
Such events are then vetoed. We are also able to reduce contributions from top production 
where one W — > jj decay occurs. 

We have employed the programs MAD GRAPH 33j and HELAS 34j to evaluate matrix 
elements for both signal and background processes. 



VI. THE DISCOVERY POTENTIAL AT THE LHC 



Based on the cuts defined above we show in Figure 3 the signal and background cross 
sections for an integrated luminosity L = 30 fb _1 and L = 300 fb _1 . The signal is shown for 
tan j3 = 10 and 50, with a common mass for scalar quarks, scalar leptons, gluino, and the 
/j, parameter from the Higgs term in the superpotential, rriq = rrig = = /i = 1 TeV. All 
tagging efficiencies and K factors discussed above are included. 

From this figure we note that the cross section of the Higgs signal with tan (3 ~ 50 can 
be much larger than that of the physics background after acceptance cuts. The Drell-Yan 
processes make the major contributions to the physics background for Higgs mass < 180 
GeV, but ti contributions become dominant for higher masses. The W + 2j contribution is 
very effectively controlled by the b tagging requirement. 

We define the signal to be observable if the lower limit on the signal plus background is 
larger than the corresponding upper limit on the background [35, 51, namely, 



L(a s + a b ) - NJL(a s + a b ) > La h + NJLa, 



which corresponds to 



(5) 



<7. s > 



AT2 



1 + 2JLa b /N 



(6) 



Here L is the integrated luminosity, a s is the cross section of the Higgs signal, and a b is 
the background cross section. Both cross sections are taken to be within a bin of width 
±AM TT centered at M^. In this convention, N = 2.5 corresponds to a 5a signal. We take 
the integrated luminosity L to be 30 fb _1 and 300 fb _1 (25| . 

For tan/? > 10, Ma and Mh are almost degenerate when Ma ^ 125 GeV, while Ma 
and Mh are very close to each other for Ma ^ 125 GeV. Therefore, when computing the 
discovery reach, we add the cross sections of the A and the h° for Ma < 125 GeV and those 
of the A and the H° for M A > 125 GeV [37 . 



Figure 4 shows the 5a discovery contours for the MSSM Higgs bosons where the discovery 
region is the part of the parameter space above the curves. We have chosen Msusy = mg = 
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(a) L = 30 fb" 1 



(b) L = 300 fb" 1 
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FIG. 3: The signal cross section at the LHC for luminosity L = 30 fb" 1 and 300 fb -1 , as a function 
of Ma, and tan/3 = 10,50. Also shown are the background cross sections in the mass window of 
Ma ± AM TT . We have applied K factors, acceptance cuts, and efficiencies of b,r tagging and 
mistagging. 



m- g 



mi = fi = 1 TeV. If Msusy is smaller, the discovery region of A , H° — > t + t~ 
will be slightly reduced for Ma ~ 250 GeV, because the Higgs bosons can decay into SUSY 
particles [38|] and the branching fraction of <p ~^ t + t~ is suppressed. For Ma ^ 125 GeV, the 
discovery region of H° — > t + t~ is slightly enlarged for a smaller Msusy, but the observable 
region of h° — > r + r" is slightly reduced because the lighter top squarks make the H° and 
the h° lighter; also the H°bb coupling is enhanced while the h°bb coupling is reduced [37 . 

We find that the discovery contour even dips below tan/5 = 10 for 100 GeV < Ma < 
300 — 400 GeV depending on luminosity. Below tan/5 = 10 our approximation of mass 
degeneracy of MSSM Higgs bosons breaks down; therefore we include only one Higgs boson 
(A ) in our calculations. 



VII. THE DISCOVERY POTENTIAL AT THE FERMILAB TEVATRON 



To study the discovery potential of this channel at the Fermilab Tevatron Run II, we 
require 

• one b quark with pr(b) > 15 GeV, \r](b)\ < 2.5 and a tagging efficiency e b = 60%, 

• one isolated lepton with pr(£) > 10 GeV and \r](£)\ < 2.0, 
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250 500 750 1000 

M A (GeV) 

FIG. 4: The 5a discovery contours at the LHC for an integrated luminosity (L) of 30 fb" 1 , and 
300 fb" 1 in the (M A , tan /3) plane. The signal includes <j) = A and h° for M A < 125 GeV, and 
(f) = A and H° for Ma > 125 GeV except that, for tan (5 < 10, (j) = A only. The discovery region 
is the part of the parameter space above the contours. 

• one jet with pr(j) > 15 GeV and \rj(j)\ < 2.5 for the tau jet, and a tagging efficiency 
of 38%, 

• the transverse missing energy ( Ifi T ) should be greater than 20 GeV, 

• the transverse mass of the lepton and missing transverse energy, Mt(1, $t)-> should 
be less than 30 GeV, 

• the transverse angular separation of the lepton and tau jet, 4>(£,j), should be less than 
170°, 

• the energy fractions for the lepton and the tau jet should be between and 1, (0 < 
X£,x h < 1), and 

• the invariant mass of the reconstructed tau pairs should be within the mass window 
of the Higg mass with AM TT = 0.1 5M 

In Figure 5 we show the signal and background cross sections for the Fermilab Tevatron. 
The signal is shown for tan/3 = 10 and 50, with a common mass for scalar quarks, scalar 
leptons and the gluino nig = nig = mi — \i — 1 TeV. All tagging efficiencies and K factors 
discussed above are included. 

From this figure we note that while hrr and ti make major contributions to the physics 
background at the LHC, brr and Wjj become the dominant background at the Tevatron 



9 



10 



100 



200 



300 



400 



M.(GeV) 



FIG. 5: The signal cross section at the Fermilab Tevatron Run II as a function of Ma, and 
tan/3 = 10, 50. Also shown are the background cross sections in the mass window of Ma ± AM rT . 
We have applied K factors, acceptance cuts, and efficiencies of 6, r tagging and mistagging. 

for Ma < 400 GeV. The cross section of the Higgs signal with tan /3 ~ 50 can be much larger 
than that of the physics background after acceptance cuts. 

Figure 6 shows the 5<r discovery contours for the MSSM Higgs bosons where the discovery 
region is the part of the parameter space above the curves. 

We find that the discovery contour for the Tevatron Run II can be slightly below tan (3 = 
30 with an integrated luminosity of 2 fb _1 and below tan/3 = 20 with L ~ 8 fb -1 for 
Ma < 150. For tan/5 ~ 50, the Tevatron Run II will be able to discovery the Higgs bosons 
up to Ma ~ 200 GeV with L = 2 fb" 1 and up to M A ~ 250 GeV with L ~ 8 fb -1 . 
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FIG. 6: The 5a discovery contours at the Fermilab Tevatron Run II for an integrated luminosity 
(L) of 2 fb" 1 , 8 fb" 1 , 16 fb" 1 in the Ma versus tan/3 plane. The signal includes cf> = A and h° 
for M A < 125 GeV, and <f> = A and H° for M A > 125 GeV except for tan < 10 where cf> = A 
only. The discovery region is the part of the parameter space above the contours. 

VIII. CONCLUSIONS 

The tau pair decay mode is a promising channel for the discovery of the neutral Higgs 
bosons in the minimal supersymmetric model at the LHC. The A and the H° should be 
observable in a large region of parameter space with tan/5 > 10. In particular, Fig. 4 shows 
that the associated final state of bcf) — ^ br + r~ could discover the A and the H° at the LHC 
with an integrated luminosity of 30 fb -1 if Ma ~ 800 GeV. At a higher luminosity of 300 
fb -1 , the discovery region in Ma is easily expanded up to Ma = 1 TeV for tan (3 ~ 50. 

In Figure 7, we compare the LHC discove ry p otential of bcfP production for the muon 
pair discovery channel, as determined in Ref.[17j|. and the tau pair discovery channel, for 
an integrated luminosity of 30 fb" 1 . It is clear that the tau pair channel can be discovered 
in a larger region of the parameter space. However, the muon pair channel can also be 
observable in a significantly large region. In addition, the muon pair channel will provide a 
good opportunity to precisely reconstruct the masses for MSSM Higgs bosons. The discovery 
of the associated final states of bcf) — > 6t + t~ and bcf) — > bfi + fi~ will provide information about 
the Yukawa couplings of bbcp and an opportunity to measure tan (3. The discovery of both 
cf> — > t + t~ and cp —>■ will allow us to study the Higgs Yukawa couplings with the 

leptons. 

We find that the discovery contour for the Tevatron Run II can be slightly below tan (3 = 
30 with an integrated luminosity (L) of 2 fb -1 and below tan/5 = 20 with L ~ 8 fb -1 for 
Ma < 150. For tan/5 ~ 50, the Tevatron Run II will be able to discover the Higgs bosons of 
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MSSM up to M A ~ 200 GeV with L 
The inclusive tau pair channel (</>° 



2 fry 



T T 



the CMS 26 



1 and up to M A ~ 250 GeV with L ~ 8 fb _1 . 
f ) has been studied by the ATLAS 19|, 25] and 



collaborations with realistic simulations. Both collaborations have confirmed 
that this channel will offer great promise at the LHC. Our results for b<p° — > br~r + are 
consistent with those given in these references and also with the results of Ref. 3^] for the 
Fermilab Tevatron. 
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